Formation of ultrashort pulses with sub-Poissonian 
photon statistics 
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A simple method for the production of ultrashort light pulses (USPs) with suppressed photon fluctu- 
ations is considered. The method is based on self-phase modulation (SPM) of an USP in a nonlinear 
medium (optical fibre) and subsequent transmission of pulse through a dispersive optical element. 
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Ultrashort light pulses (USPs) continue to draw the 
attention of investigators today. The state of the art in 
this area of laser physics in the late 1980s was set forth 
by S. A. Akhmanov et al. in Refs. 0,0- In the past 
decade, considerable progress has been made in quan- 
tum optics in the generation of nonclassical (the so called 
squeezed) light fields. The formation and application of 
USPs in a nonclassical state makes possible to combine 
in experiments a high time resolution with a low level of 
fluctuations. 

In principle, one can obtain pulsed light fields in a 
nonclassical state by using the same nonlinear optical in- 
teraction as those used in the case of continuous fields || , 
Parametric amplification is a technique that is most 
extensively used for this purpose nowadays. In the case 
of degenerate three-frequency parametric amplification, 
quadrature-squeezed light is produced. However, this 
light is found to have super-Poissonian statistics directly 
at the output of the amplifier, and one needs interfer- 
ometers to transform it to get sub-Poissonian statistics. 
One can obtain light with sub-Poissonian photon statis- 
tics with the aid of nonlinear interferometric devices in 
the presence of self- phase modulation (see, e.g., Rcf. ||). 
Note that the self-phase modulation itself is not accom- 
panied by a change in photon statistics. 

In a recent paper, we studied SPM of a light pulse and 
its subsequent propagation in a dispersive linear medium 
(or the passage through optical compressors). In this 
case, pulses with sub-Poissonian photon statistics can be 
formed. We managed to make an accurate calculation 
of the process under consideration owing to the consis- 
tent quantum theory developed by us for the self-action of 
light pulses in a medium with inertial nonlinearity || , || . 
The method of formation of pulses with sub-Poissonian 
statistics proposed here is simple for experimental real- 
ization and stable against external random effects (tech- 
nical fluctuations). Note that the theory developed in 
Refs. ||, j(| for the self-effect of USPs takes into account 
a relaxation time of a nonlinearity that determines the re- 
gion of the spectrum of quantum fluctuations which are 



of substantial importance in the formation of squeezed 
light and does not limit the amplitude of quantum fluc- 
tuation. 

When analysed from the quantum point of view, SPM 
of USPs is described by the expression 0j , M . 



A(t,l) = e d ^A (t), 



(1) 



and the Hermite conjugate operator A + {t, I). Here 
A + {t, I) (A(t, I)) is the photon creation (annihilation) op- 
erator for the given section i at a given moment of time 
t (the output of a nonlinear medium is specified by the 
section x = I, and the input is specified by the section 
x = 0, A(t,x=0) = A (t)), 6(t) = ijq[no(t)], and 7 is the 
coefficient related to nonlinear properties of a medium 
and proportional to the length / in it. The expression 
7<7[rio(t)] characterizes the nonlinear phase incursion, 



q[Mt)} = / ff(|ti|)n (t-ti)dti, 



(2) 



and no{t) = Aq (t)Ao(t) is the operator of the "density" 
of the number of photons at the input of a nonlinear 
medium. The function H(t) takes into account a finite 
response time of the nonlinearity of a medium. In view of 
the causality principle, H(t) ^ for t > and H(t) = 
for t < 0. In expression (|]) and (0), t = t — x/u is the 
time in the moving coordinate system, t is the current 
time, and u is the pulse velocity in a medium. 

In the case of SPM (see (□)) discussed here, the 
quadrature-squeezed light || , || is formed, and the pho- 
ton statistics in a medium is not changed because the op- 
erator of the number of photons n(t,l) = ho{t) remains 
unchanged. In what follows, we show that the propa- 
gation of a pulse (see ([l])) through a dispersive optical 
device is able to transform its initial Poissonian statistics 
into sub- or super-Poissonian statistics. 

For example, let us consider the propagation of a pulse 
in a dispersive linear medium. At the output of the 
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medium, we have the following expression for the pho- 
ton annihilation operator H: 



(3) 



B(t,z)= / G(t-t 1 ,z)A(t 1 ,l)dt 1 . 



It is known that the operators A(t, I) and B(t, z) 
are bound to satisfy commutation relation of the 
[C(ti), C + (t 2 )] = S(ti — t?) type. In view of relation 
(H) between the operators under consideration, we find 
the condition imposed on the Green function G{t) of the 
dispersive clement: 



G{t x - t, z)G*(t 2 - t, z)dt = S(t 2 - h). 



(4) 



Let us introduce the operator of the number of photons 
over the measurement time T, 



t+T/2 

N T (t,z)= J N(tx,z)dh, 
t-r/2 

and define the Mandel parameter 

e(t,z) 



(5) 



Q(t,z) 



(N T (t,z)) 



(6) 



where 



N(t,z) = B+(t,z)B(t,z); 

e(t, z) = (N*(t, z)) - (N T (t, z)) 2 - (N T (t, z)). (7) 

The parameter Q(t, z) characterizes the difference be- 
tween photon statistics and Poissonian statistics (for the 
latter Q(t,z) = 0). The angle brackets in (0) denote 
averaging over the initial quantum state of the pulse. 

One can calculate expression (R) by using the algebra 
of time dependent Bose operators || , B . In this case we 
have the following relations: 



O o(t) 



Ncxp 



ho{t — txjdtx 



(8) 



(9) 



where TL{t) = i-fh(t), h{t) = H(\t\) and N is the operator 
of normal ordering. 

In the case of a coherent initial pulse, a small nonlinear 
phase incursion per photon, and the pulse duration r p is 
considerably grater than the times T and r r (r r is the 
nonlinearity relaxation time). After some cumbersome 
calculations, taking into account expression (Q) and (j^) 
one gets 



:(t,z) = 2jT z Im[lf(t,z)r 2 (t,z)} 



where 



h(t,z)= / G(t-ti,z)a(ti)e l,/ ' (tl) ^i; 



(10) 



(11) 



h{t,z)= J J G{t-t 1 ,z)G(t-t 2 ,z)a{t 1 )a{t 2 ) 

— oo — oo 

x h(h - t 2 )e l W^+^^dtxdh; (12) 

Here a(t) is the eigenvalue of the operator A (t); ifj(t) = 
2jno(t) is the nonlinear phase incursion, n-o(i) = |«(t) | 2 
is the average density of the initial number of photons, 
and (N T (t,z)} =T\h(t,z)\ 2 . 

Let us assume that the Green function G(t) has the 
form 



G(t,z) = {-i2irk 2 z)- 1/2 



exp 



if* 
2k 2 z 



(13) 



where z is the distance traveled in the dispersive linear 
medium. We have the coefficients k 2 > for normal dis- 
persion of the group velocity, and k 2 < for anomalous 
dispersion. Expression (|l3|) corresponds to the case when 
dispersive properties of a medium are taken into account 
in the second approximation of the dispersion theory [pj . 

The dispersion effect is known [Q to have a spatial ana- 
logue in the form of diffraction spreading of a wave beam. 
The quantum description of diffraction in nonlinear opti- 
cal processes was considered first by S. A. Akhmanov et 
at in Ref. ||] (see also Ref. ||), where it was shown that 
the mixing of different angular components of a para- 
metrical amplified beam is able to cause additional non- 
classical effects. One can expect similar effects in the 
case under consideration. The self-interaction effect of a 
pulse given by expression (|l|) is accompanied by changes 
in its phase (and therefore frequency) over time. In the 
course of its passage through dispersive linear elements, 
a frequency-modulated pulse of this type is compressed 
or stretched j2|. It is precisely this effect that is able to 
change its photon statistics. 

Let us determine the parameter Q(t,z). To simplify 
calculations, we take the nonlinearity response function 
in the form h(t) = (l/r r ) exp (— i 2 /r 2 ). Let a pulse be 
initially of the Gaussian form n (t) = n exp (— t 2 /r 2 ) . 
In the paraxial approximation, i.e., in the case where the 
phase tjj(t) is replaced with ^(l — t 2 /r 2 ), (ip = 2^n ), 



we obtain for r p > r r 



(N T (t,z))=n TV- 1 (z)ex.p 



V 2 (z)i 



(14) 



<3(o, z) 



-tpo 



T 



sin< arctan 



<P( Z ) 
w(z) 



2 



FIG. 1. Dependence of the parameter Q(0,z) [Z] in the 
case of compression of a USP with SPM on the nonlinear 
phase incursion ipo [X] and the dispersion phase <p(z) [Y] for 
r p /T = 8. 



1 / 2ip d (z)w(z) 

+ 2 arCtan U^)^W"^(z) 

I [vo\z) - 2<p 2 (z)m 2 (z) + lf \ (15) 



where 



V 2 (z) = m 2 (z) + ip 2 (z); w(z) = 1 - s^ ip{z); 

D d |fc 2 | \ki\ 

D and e? are the characteristic dispersion lengths, s — 1 
for k2 < 0, and s = — 1 for fe 2 > 0. Let us restrict our 
consideration to the analysis of the parameter Q(t, z) for 
t = 0, because expression (^5j) for an arbitrary moment 
of time t is rather cumbersome. 

From ( |l4| ) and (|lq ) it follows that one can have 
z) < 0, i.e., obtain a pulse with sub-Poissonian pho- 
ton statistics through both compression (s = 1) and 
spreading (s = — 1) of a pulse. Of particular interest from 
the practical point of view is the case of anomalous group 
velocity dispersion (s = 1) in which the compression of 
a phase-modulated pulse (see (H)) takes place. The de- 



pendence of the Mandel parameter and the average num- 
ber of photons in this case are presented in Fig]l]-Fig.H 
and Fig.||. One can see that the suppression of fluctu- 
ation of the number of photons becomes noticeable for 
the nonlinear phase ipo > 1. For a given nonlinear phase 
■00 ; there is a certain dispersive phase incursion for which 
the suppression of the number of photons is the greatest. 
According to (|l4|), this is obtained for ip pt{z) = l/i/V ln 
this case, the Mandel parameter has the minimum value 
Qmin = — [1~/t p ]iI> 2 . Note that a pulse has the mini- 
mal duration for tp(z) = 0o/(l + ipo)- For Vo S> 1 both 
extremities are almost coincident. 



FIG. 2. Dependence of the parameter Q(0,z) [Y] in the 
case of compression of a USP with SPM on the dispersion 
phase <p(z) [X] for tpo = 4 and r v /T = 8. 



As follows from calculations, the photon statistics for 
the measurement time T ^> r p is of the Poissonian type. 
A situation identical in many aspects to the one con- 
sidered here was observed in spectral measurements of 
fluctuations of femtosecond pulses transmitted through 
an optical fibre in the region of normal group-velocity 
dispersion. An increase in the spectral bandwidth of a 
filter caused an increase of the level of fluctuations || . 

Let us present a numerical example. Consider the in- 
teraction of a pulse at a wavelength of 1 fim in an optical 
fibre, and let its initial duration and maximum inten- 
sity be 2r p = 10 ps and 10 7 W ■ cmT 2 . In this case, 
the nonlinear phase in a quartz fibre for the path length 
/ = 100 m is ipo = 3. If the resulting USP travels through 
an optical compressor with anomalous dispersion fc 2 | = 
10 -26 s 2 ■ cm" 1 , the dispersion phase (p pt(z) = 1/00 
is obtained at a distance of about 8 m. This exam- 
ple illustrates the feasibility of formation of USPs with 
sub-Poissonian photon statistics by the method proposed 
above. 

The obtained result shows that it is possible in princi- 
ple to form a pulse with sub-Poissonian photon statistics 
in the case where a high-intensity pulse undergoes SPM 
during its travel through a nonlinear medium and sub- 
sequently passes a dispersive optical device (an optical 
compressor or an optical fibre) . In view of the fact that 
these processes occur in succession, one can choose an 
optimum scheme or optimum conditions at each trans- 
formation stage. In the specific case of formation of low- 
intensity USPs in a nonclassical state, one can atten- 
uate high-intensity pulses at the output of a nonlinear 
medium. It is natural that this is accompanied by a par- 
tial loss of nonclassical properties. It should be noted 
that in the case of formation of optical solitons, SPM 
and the dispersion effect take place simultaneously . 

Authors are grateful to V. A. Vysloukh for fruitful discus- 
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FIG. 3. Variation of the shape of a USP with SPM dur- 
ing compression. The calculation are made in the paraxial 
approximation. ([Z] = (N r (t,z)), [X]=t/r p , [Y]= <p{z)). 



sions. This work was supported in part by the Fundamental 
Metrology Program of the State Committee on Science and 
Technology. 



[1] S. A. Akhmanov, V. A. Vysloukh, A. S. Chirkin, Usp. Fiz. 
Nauk., 149, No.450 (1986). 

[2] S. A. Akhmanov, V. A. Vysloukh, A. S. Chirkin, Optica 
Femtosekundnyh Lazernykh Impulsov (Optics of Femtosec- 
ond Laser Pulses) , Moscow, Nauka, 1988 (New York: AIP, 
1992). 

[3] S. A. Akhmanov, A. V. Belinski, A. S. Chirkin, in 
Novye Fizicheskie Prvntsipy Opticeskoi Obrabotki Infor- 
matsii (New Optical Methods of Optical Information 
Processing). (Eds S. A. Akhmanov, M. A. Vorontsov) 
(Moscow, Nauka 1990) pgs. 83-194. 

[4] D. F. Walls, G. J. Milburn, Quantum optics, (Berlin: 
Springer, 1995). 

[5] F. Popescu, A. S. Chirkin, Pisma Zh. Eksp. Teor. Fiz., 
69, 482 (1992), [JETP Lett., 69, 516 (1999)]. 

[6] F. Popescu, A. S. Chirkin, quant-ph/ 0003028/ 8.03.2000/ 
Submitted to Phys. Rev. A. 

[7] S. A. Akhmanov, Yu. E. Dyakov, A. S. Chirkin, Vvede- 
nie v Statisticheskuyu Radiofiziku i Opticu (Introduction 
to Statistical Radiophysics and Optics) (Moscow: Nauka 
1981). 

[8] S. A. Akhmanov, A. V. Belinski, A. S. Chirkin, Kbanto- 
vaya Electron. (Moscow) 15, 873 (1998), [Quantum Elec- 
tron. 18, 560 (1988)]. 

[9] F. Konig, S. Spatter, I. L. Shumay, A. Sizmann, Th. 
Fauster, G. Leuchs, J. Mod. Opt, 14, 2425 (1998). 



